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I.  INTRODUCTION. 

Upper  bounds  for  multivariate  normal  probability  areas  outside  of 
rectangles  centered  at  the  origin  are  of  interest  due  to  their  applica¬ 
tions  in  producing  conservative  unbiased  simultaneous  confidence 
intervals  and  hypothesis  tests.  Unfortunately,  it  is  often  not  compu¬ 
tationally  practical  to  integrate  multivariate  normal  area  over 
dimensions  n  >  A  and  obtain  these  probabilities  exactly.  The  method 
commonly  used  to  determine  upper  bounds  fot  normal  probabilities 
of  these  regions  is  based  on  the  conservative  assumption  of  indepen¬ 
dence  as  given  in  the  following  theorem  suggested  by  Dunn  (1958)  and 
proved  by  Sidak  (1967). 

Theorem  1-1.  Let  x  '•  N(0,~)  and  let  A.  be  the  event  x  .  d  (-C  .  , c  .  )  '  . 

-  -  i  ill 

Then 

n  n  n  n 

Pr{  -  A? '  <  r.  PriA?}  ->  Pr{  A.}  >  1  -  7  Pr{A?i  . 

1  '  -_i  1  .  ,  i  —  .  -i  i 

i=l  i=l  i=l  i=l 


Since  1967,  three  different  approaches  have  been  taken  which  enable 

n 

one  to  obtain  lower  upper  bounds  for  the  Pr{  A.)  than  those  given  bv 

i=l  1 

Theorem  1-1.  These  approaches  require  that  £  be  known  and  integration 
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over  joint  m-variate  densities  be  possible  for  some  m:  1 


<  m  <  n 


to  obtain  Pr-.A.  A.  u***u  A.  }.  These  three  approaches  are  described 
1  12  1m 

in  chapter  II,  and  compared  in  chapter  III.  They  are  combined  together 

n 

in  chapter  IV  to  produce  lower  upper  bounds  for  Pr{  ij  A.}  than  are 

i=l  1 

obtained  from  using  any  one  approach  alone.  In  chapter  V  examples  of 

combining  these  methods  are  given  and  in  chapter  VI  selected  tables 

n 

of  upper  bounds  for  Pr{  •,  i  A.}  given  different  combinations  of  the 

i=l 

three  approaches  are  shown.  Chapter  VII  has  a  descriptive  summary  of 
the  terminology  used  in  this  report  and  can  be  u=d  for  quick  reference. 
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II.  NEW  APPROACHES  FOR  UPPER  BOUNDS. 

Approach  I.  Intersection  Subtraction 

This  theorem  was  first  proven  for  m  =  2  by  Hunter  (197b) 
and  later  extended  to  m  >  2  by  Hoover  (1987).  The  theorem  is 

Theorem  A-I.  Let  x  be  a  n v  1  vector  distributed  N(0,I),  y 
be  any  permutation  r  ( • )  of  the  elements  of  x,  and  A^  he  the 
event  v.  t  (-c  ,.,,c  , ..).  Let  1  <  m  <  n  and  let  s.  be  a  set 
of  size  (m-1 a  i-1)  of  integers  l,2,...,i-l.  Then 


Pr {  x . i (-c .,c.)} 

i=l  1  1  1 


Pr-.A.  )  +  Pr •' A .  ;  -  Pr'A. 
1  l  i 


i=2 


-i£si 


A.), 


A  . 


=  Pr .A,  '  +  Pr;A.  r  ( 

1  i=2  1  jcs.J 


Proof . 


Pr{  '  x . t  (-c  .  ,c . ) 1  =  Pr{  ('  y . t  (-c_  , . .  ,c_  ) 

i-1  111  i=!  1  <x) 


n 


n 


i-1 


=  Pr;  (A.);  =  Pr-.A  ;  +  PrlA.  “  [  A. 
i=l  1  i=2  1  j=l  J 


n  n 

Pr'A  }  +  7  Pr { A  .  I  i;  A.]C>  =  Pr(A  } +  )  Pr ( A , '-Pr  \  A .  [  :  ,  A.] 

1  i=2  1  jts.  1  i=2  1  1  jcs  J 

l  i 


Remarks . 

1)  To  implement  this  theorem  requires  being  able  to  determine 

Pr'A.  [  i  A.]C(,  which  requires  being  able  to  integrate  over  m 

j^s.  J 

dimensional  marginal  distributions. 
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2)  For  m  =  2  the  above  theorem  is  equivalent  to  Hunter's 
Theorem  which  is:  Let  x  be  0  N(0,I).  Let  I  be  any  maximal 
tree  connecting  the  x_.  and  forming  no  circuit  and  e^  be  edges 
on  this  tree  then 

n 

Pr{  \)  xi^(-c±,c  .) }  <_ 
i=l 


'  Prfx . t  (-c . ,c  . ) } 

l  i’  l 


\  Pr(x . t  (-c . ,c . )  0 x .^(-c . ,c . ) 3 . 

e  eT  1  11  J 

ij 


3)  This  theorem  still  leaves  unsolved  how  to  choose  the  "best" 


permutation  -(•)  of  x^  and  sets  s^  to  give  the  lowest  upper 
n 


bound  for  Pr{  x  ,i(-c .  ,c .)  .  At  least  one  "best"  permutation  set 
i=l  1  1  1 

s^  exists  since  there  are  a  finite  number  of  permutations  and  sets 
s^.  For  the  case  m  =  2  this  is  equivalent  to  finding  the  best  tree 
T  for  the  theorem  mentioned  in  remark  2.  This  is  possible  to  <->o; 
the  procedure  is  given  by  Theorem  AP-1  in  the  Appendix. 

4)  This  approach  is  distribution  free.  The  theorem  is  true  for 
all  unions  of  any  collection  of  events  A^ :  i  =  l,...,n.  The  fact  that 
A^  was  (i) ,c- (i)  ^  where  y  ^  "  Nv,0,r)  and  “(•)  is  a  permu¬ 

tation  of  l,...,n;  was  not  used  in  the  proof.  The  next  two  approaches 
will  be  less  general. 


Corollary  A-IC.  Group  Subtraction  Overlap. 


Let 

si 

and 

s 

^  be  sets  of  integers  from  ; 

■  1,  • • • -n' 

Let 

h 

be 

the 

upper  bound  for 

Pr {  ,  A. } . 
.  '  l 

1ES1 

Let 

*  7 

be 

the 

upper  bound  for 

Pr;  j,  A. ; 

.  l 

itS2 

and  s^. 

Let 

S3 

be 

the 

intersection  of 

the  sets  s^ 

-4 


Then 


Pr{  y  A.}  <  {.  +  L.  -  Pr {  .  ,  A.} 

.  /  .1  —  1  2  w  l 

ltCs^us,,)  ies^ 


Proof . 


pr{  u  At>  =  pr{  U  A  }  +  pr{  y  A.}  -  Pr{  (  y  A .  )<.  (  y  A.)} 
ie(s^as2)  ies^  1  ies2  ^  ies^  1  ies2  J 

1  pr{  y  A.}+Pr{y  A . }  -  Pr  {  y  A  ’ 

ies  1  ies2  J  icCs^s^ 


Pr{ 


its 


A. }  . 

i 
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Approach  II.  Conditional  Multiplicative  Approach. 

The  following  theorem  is  a  corollary  of  a  theorem  from  Glaz  and 
Johnson  (1984), 

Theorem  A-II.  Let  x  be  an  n>  1  vector  distributed  N(C,T)  where 
all  off  diagonal  elements  of  -Dl'D  are  nonnegative  for  some  D, 
where  D  is  a  diagonal  matrix  with  elements  +1.  Let  v  be  any 
permutation  *'(•)  of  the  elements  of  x  and  A^  be  the  event 
v.c(-c  ,.,,c  Let  1  <  m  n  and  let  s.  be  a  set  of  size 

(m-lAi-1)  of  integers  i,2,,,.,i-l  for  i  =  2,...,n. 

Then 


n 

Pr {  •  .  x .  t  (-c . , c . ) } 
.  ,  i  i  i 

i=l 


n 

1-2  Pr { A 
i=l 


AC 1 

1  l 

l„s^ 
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Proof . 


Part  One. 


n  n  n 

Pr{  U  x.£(-c.,c.)J  =  1  -  Pr{  n  'x.e(-c . ,c . ) }  =  1-Pr{  ^  (y . c (-c  , . .  ,c_  . . . 
.  .  1  1’  1  .  '  i  i  i  .  ,  i  tt(i)  d) 

i=1  i=l  i=l 


1  ~  i=1Pr{yiC (-C7T (i) ,C^ (i)  '>’jc('cn(j)’cn(j))’  j 


Part  Two.  By  a  theorem  in  Karlin  and  Rinott  (1980) 

(a)  jv!  has  an  MTP^  distribution  when  xiN(0,I)  and  the 

elements  of  -DI  are  nonnegative  where  D  is  any  diagonal  matrix 

with  elements  +1  . 

(b)  (-o^(i)  ,c_  (_.))  for  all  i“'  j  j  e  [  0  ,c_  ^ )  for  all  i  and 
these  are  monotone  sets  of  the  same  type. 

It  hence  follows  from  (a),  (b)  and  theorem  2.3  part  (1)  in  C.laz  and 
Johnson  (1989)  that 


Pr{yiC(-c_(i)iy.c(-c_(.),c_(.))  for  all  jcs.}  < 


Pr{yi£(-cT(i),c_(i))|yjE(  c_  (j ) ,c, (j } ) ,  j*l, . . . ,i-l: 


which  implies  that 


|-c-u)'c”(i)>:!]‘K(i)’c-(i))  for  aU  3!si'- 


i=l 


Pr{yiE(-c._(i),cr(.))  'vj:(-c_(.),c^(.)),  j  =  l - ,i-l  ' 


See  Karlin  and  Rinott 


1980  for  a  definition  of  MTP?. 

L 
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which  when  applied  to  the  result  of  part  one  implies 


Pr^Xi^-Ci.c.)}  <  l-iriPr{yic(-c_(i),cr(i))iyje(-c.(j),c;.(j)),  jts.} 


=  1  -  r.  Pr {a?  !  U  Ac} 
♦  1  1  .  1 
1=1  JCS.  J 

1 


Remarks . 


1)  To  implement  this  theorem  requires  being  able  to  determine 

Pr(A^ ;  A?}  which  requires  being  able  to  integrate  over  m 

jts.  J 

dimensional  marginal  distributions. 

2)  For  2  the  above  theorem  is  equivalent  to  the  following: 

Let  x  -  N(0,._)  where  the  off  diagonal  elements  of  -DTD  are 

nonnegative  when  D  is  some  diagonal  matrix  with  elements  +1. 

Let  T  be  an>  maximal  tree  connecting  the  x^  and  forming  no  circuit 

WLOG  let  be  a  terminal  mode  of  this  tree  and  let  all  connections 

be  directed  toward  .  Let  e..  be  the  directed  edge  connecting  x.  tc 

1  ij  .1 

e  . 

•  i  i  c 

x.  i.e.  (x.  _  x.)  and  let  A.  be  the  event  x.c(- c.,c.).  Then 

i  j  i  l  ill 


■{  l  x.tf (-c.  ,c. ) }  >  1  - ( Pr (A,  }C  H  Pr  "A? | A? } ] . 
11-  1  _ 


e.  .  T 
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3)  As  in  Theorem  A-I  for  m=2  the  best  tree  T  can  be  found  using 
the  procedure  of  theorem  AP-1  in  the  Appendix.  For  m  >_  3  it  is  harder 
to  find  tlio  "best"  -permutation  (•)  and  sets  to  get  the  lowest 
possible  upper  bounds  from  this  method. 

4)  This  theorem  is  not  distribution  free  it  depends  highly  on  the 


fact  that  !xi  is  mtp  and  that  intervals  of  the  form  [O.c.)  are  monotone 

~  2  l 

of  the  same  type. 


/ 


Approach  III.  Independent  Subunit  Conservative  Bounds. 

Theorem  A-III.  (Khatri  1970.  Proved  by  Bechner  1987).  Given  x  N(0  ') 
and  { ),  {h'2 {Wj.}  are  disjoint  collections  of  such  that 


F 

(W  }  =  ;x  ,x? 
f=l  1 


Then 


n  F 

Pr\  x.c(-c.,c.)}  >  7  Prt  ^  x.c(-c.,c.)} 

..1  i’i  -  ,  ,  .,  x  11 

1=1  f=l  Xi-WF 

which  implies 

n  F 

Pr\  '  x  .  i‘  (-c  .  ,c  . )  ■  1  -  7  Pr{  '  x.c(-c.  ,c.)  •. 

i=l  f=l  x . ek 

1  F 

Proof .  The  proof  in  Khatri' s  article  was  incorrect.  It  has  come  to  my 

attention  that  the  theorem  lias  been  subsequently  proven  by  Bill  Bechner 

at  the  University  of  Texas  Austin  but  the  proof  has  not  been  published  yet. 
Remarks . 

1)  By  itself,  this  theorem  is  not  any  superior  to  the  Dunn-Sidak 
method  but,  as  we  shall  see  later,  is  useful  when  combined  with  the  two 
previous  approaches  mentioned  in  this  chapter. 

2)  This  theorem  is  highly  nondistribution  free  and  can  be  used 
only  on  convex  symetric  probability  regions  centered  at  the  origin 

(or  compliments  of  such  regions)  involving  densities  where  f(x)  -  f(-x). 

Corollary  A-II1C.  Given  x  N(0,2)  and  {W^ (V  j  are 

disioint  collections  of  the  x.  such  that 
J  1 


F 


f=l 


tlv. 


lXX  »x2 ’ 1  * ' 


,x 
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and  v  r  is  an  upper  bound  for  Pr 

i  S  -&.1  (-c  .  ,  c  .  ) }  .  Then 
l  l  l 

icUf 

n 

F 

Prl  )  x  .  i  (-c  .  ,c  . ) } 
i=l  1  1  1 

-  i  -  n  (i-i  )  . 
f=i 

Proof . 

B\  Thm.  A 
n 

-hi  F 

P r {  ;  x.^(-c.,c.)}  < 

i=l  1  1  1 

1  _  pr  (  n  x.c(- 

f=l  itW  1 

F 

=  1  -  :  1  -  Pr;  V  x 

C-C  ,G  }  ■ 

f =1  i::Wf 

X  X  J. 

F 

l  -  "  d-:f)  • 

f=i 

bv  repeated  application 
of  the  distributive 

multiplication  inequality 
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4*. 


III.  COMPARING  THE  METHODS. 

Thee  eras  A-I  and  A-II  have  been  presented  in  a  way  to  enable  us 
to  use  m  dimensional  marginals  for  m  2.  To  keep  the  statements 
and  proofs  of  the  theorems  simple  for  the  remainder  of  this  technical 
report  m  =  2  will  be  assumed  whenever  theorems  A-I  and  A-III  are  used. 

All  of  the  results  in  this  chapter  could  be  extended  to  comparable 
results  for  m  _>  3. 

Sine'-  there  are  several  different  methods  to  use, the  question 
arises  as  to  under  which  conditions  is  one  method  better  than  the 
others  (i.e.,  gives  closer  bounds).  This  question  is  explored  further 
in  the  present  chapter. 

Thereon  C-l.  Let  x  and  A^  be  as  described  in  Theorems  A-I,  A-Il 
and  A-III.  l.et  n  =  2  then  (i)  the  intersection  subtraction  method 
(Theorem  A-I)  upper  bound  is  exactly  the  same  as  the  conditional  multi¬ 
plicative  method  (Theorem  A-II)  upper  bound  and  is  exact.  (ii)  When 
the  correlation  is  non  zero,  the  Dunn-Sidak  upper  bound  result  is 
inferior  to  the  upper  bound  of  theorems  A-I  and  A-II. 

Proof .  (i)  The  A-I  upper  bound  for  Pr;A^  is  Pr;.A  i+Pr.A,'- 

Pr-A  'A,,  which  equals  Pr  -  A  A0'-  by  elementary  axioms  of  probability. 

The  A-li  upper  hound  for  Pr{A^.A  ■  is 

Pr!A2  Ai’  ’  1-Pri#rA21  ■  PriM  V 

(ii)  Theorem  1  in  Sidak  (1967)  and  its  proof  show  that  Pr{ A, ■  is 

a  monoton i ca 1 ly  decreasing  ! unction  of  the  absolute  value  of  the  correla- 

n  "> 

tion  coefficient  for  fixed  variances  cT  and  where  again  A.  = Pr  x.  i  (-c.,c.l 

1  2  'i  til 
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Theorem  C-2.  Let  x  and  A^,  i=l,...,n  be  as  described  in 
theorems  A-I,  A-II  and  A-III.  For  n  >  2,  assuming  m  =  2  and  that 
the  conditions  are  met  for  the  conditional  multiplicative  (A-II) 
upper  bound,  then  for  any  upper  bound  given  by  either  overlap  subtrac 
tion  (A-I)  or  by  Dunn-Sidak  (1-1)  there  exists  an  equal  (or  superior 
if  I  is  not  diagonal)  upper  bound  given  by  (A-II). 

Proof .  First.  For  a  given  directed  tree  T  with  edges  e„  and  the 
root  at  x ^  as  described  in  remark  2  of  theorem  (A-II), 
the  (A-II)  upper  bound  is 


1 


[Pr(AC)  V.  (AC!AC)] 

1  e .  . ;  T  1  J 
i  j 


1  -  Pr(AC) 
1 


Pr  (A^ c  A?) 

e . .  cT  Pr(AC) 
ij  1 


Pr(AC)Pr(AC)  n 

__  1  -  Pr(A  )C  T - - - =  1  -  V.  P r  ( A . )  ° 

1  e .  .  T  Pr(AC)  i=l  1 

by  theorem  1  ^  ^ 

1  in  Sidak 

(The  inequality  is  strict  if  T.  is  not  diagonal) 


=  I'unn-Sidak  upper  bound  . 


Second.  For  a  given  directed  tree  T  with  edges  e..  and  the  root  at 
x  as  described  in  remark  2  of  theorems  (A-II)  and  (A-I).  Let  T,  be 
the  tree  connecting  the  first  two  elements  in  T  and  T.  be  the  sub 

l 

tree  of  I  connecting  the  first  i  elements  in  1  i=3,...,n.  Let 


(A-l)  .  and  (A-II).  be  the  A-I  and  A-II  upper  bounds  for  Pr  '  (A.) 

j  ;  Ti  1 

Then  bv  theorem  C-I  v;o  have 


(A-I)2  =  (A-II),  . 
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Now 


Thus  (A- I ) 3  ■  (A-I i ) 2  . 

following  this  proof  inductively  expressing  (A-I)^  and  ( A- 1 1 ) ^  in  terms  of 

(A-I),.  and  (A — 1 I ) 0  and  x.  the  fourth  element  of  tree  T,  we  will  get  (A-I).  - 

5  i  A  h  '■  <4 

(A— II) .  etc.  until  obtaining  (A-I)  >  (A-I I)  for  anv  n  ^  2.  But  T  =T  and  thus 

4  n  n  n 

the  upper  bound  from  using  A-I  on  T  is  higher  than  the  upper  bound 
from  using  A-I I  on  I  which  completes  the  proof. 
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theorem  C-3.  Let  T  be  as  described  in  remark  2  of  theorems  A-II 

2 

and  A-III  and  let  o.  and  c.  be  constant.  If  max'P.. 1  is  bounded 

l  l  T  ‘  ij 

by  some  value  less  than  one  then  as  n  -*■  <*,  the  Sidak-Dunn  assumption 
of  independence  eventually  produces  a  better  upper  bound  for 


Pr {  O  A. 
i-1  1 


than  does  the  intersection  subtraction  method  (A-I)  and  in  fact  for  n  large 
enough,  (A-I)  will  produce  an  upper  bound  for 


Pr{  ' A.  ) 
1=1  1 


which  is  larger  than  1. 


Proof .  The  intersection  subtraction  method  upper  bound  is 


p  r  A  •  -  Pr  A .  -  A.  '  = 

i  -  1  J 

e  .  .  ;  i 
1J 

Pr;A  •  +  ;  [Pr{A  >-Pr{A.  A  ;i]  > 

e .  ~  cT  1  J  f 

1J  *  ‘ 

Bv  Theorem  1  of  Sidak  since  P...  is 

ij 

bounded  below  1,  Pr ■ A . ; -Pr (A . ' A . )  > 

i  i  J  “ 


n > 0  for  all  i  and  j.  See  footnote  below 


i  [  r:  ]  =  Pr  ■ A  .  ;+  (n-1)  r;  ■>  1  s  l-~  (1-Pr  ’A  .  0  =  independence 
e  cl  fc  1  upper  bound 

ii  \  r _ •, _ _ u 


LJ  \ 

for  some 
>  >  o 


for  n  large  enough 


Since  c.  and  c.  are  constant  and  max  P..1  is  bounded  bv  some  value 
ii  T  '  1 1 

less  than  one,  it  follows  that  Pr'A.  ■  =  k, ,  for  k,  some  constant  and  bv 

i'll 

Theorem  1  in  Sidak  that  Pr-'A.  'A.}  <  k„  for  k„  k.  .  Therefore, 

i  J  -  2  2  1 

[Pr {A .  }-Pr{A .  ■-•A . }  ]  >  n  for  some  constant  n  >  0. 
i  i  J  ~ 
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Conclusion. 


Theorem  C-2  shows  that  whenever  it  applies,  the  conditional  multi¬ 
plicative  (Theorem  A-II)  always  gives  ti  'hter  upper  bounds  than  do 
either  the  intersection  subtraction  (Theorem  A-I)  or  the  Sidak-Dunn 
(Thereom  1-1)  methods.  Hence  (A-II)  should  be  used  in  preference 
to  (A-I)  and  (1-1)  to  obtain  upper  bounds  whenever  the  data  is  jointly 
MTP2.  Theorem  C-l  implies  that  for  small  n,  the  intersection  sub¬ 
traction  method  (A-I)  gives  tighter  upper  bounds  than  does  the  Dunn- 
Sidak  method  (1-1).  Theorem  C-3  implies  that  as  n  gets  larger,  the 
Dunn-Sidak  ( I— 1 )  upper  bounds  tend  to  become  superior  to  those  of  the 


intersection  subtraction  method  (A-I)  for  cases  with  the  same  c^. 


and  IP  .1  bounded  away  from  one.  The  Dunn-Sidak  method  has  been 
compared  with  the  intersection  subtraction  method  in  more  detail 


(see  Hoover  (1986)). 


IV.  COMBINING  THE  APPROACHES. 


The  three  methods  A-I,  A-II  and  A-III  can  be  combined  together 

n 

to  give  lower  upper  bounds  to  Pr{  y  x.£(-c.,c.)}  than  could  be 

i=l  1  11 

obtained  by  using  any  of  these  methods  ,  or  the  Dunn-Sidak  Theorem  alone  . 
In  fact,  whenever  E  is  nondiagonal,  there  will  exist  some  combina¬ 
tion  of  A-I,  A-II  and  A-III  which  gives  a  lower  upper  bound  to 
n 


Pr{  ii  x.£(-c.,c.)}  then  does  the  currently  used  Dunn-Sidak  method. 
i=l  1  1  1 

It  would  be  nice  to  be  able  to,  for  a  given  situation,  determine  the 

best  combination  of  A-I,  A-II  and  A-III  in  terms  of  giving  the  lowest 

n 

upper  bounds  for  Pr  {  ^  x . t  (-c . , c . ) } .  Unfortunately  in  most  cases 

i=l  ill 

this  would  be  computationally  impractical  to  do.  Nevertheless  it 
is  possible  to  develop  procedures  which  combine  the  methods  in  ways 
that  produce  good  result^. 


Example  M-one. 

At  this  point,  a  description  and  example  of  how  the  three  methods 
can  be  used  together  should  be  given.  Let  us  assume  x  is  an  81 
vector  with  a  N(0,T)  distribution  and  that  a  directed  tree  with  the 
following  links  has  been  set  up 


x 


4 


X2  "  X3  "  X7  x8 


Ke  now  define  three  types  of  links 

S  -  Intersection  subtraction  (A-I) 

M  -  Conditional  multiplicative  (A-II) 

I  -  Conservative  independent  subunit  (A-III) 
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These  links  correspond  to  which  method  is  being  used.  We  can  now  assign 
these  link  types  to  such  links  on  the  tree  for  instance 


4 

I  S 

SIMS 
X1  X2  *"  x3  ■"  x7  ■“  x8 
•f  M 

X5 
t  M 


The  first  step  we  wish  to  take  is  to  identify  the  continuous  sets  of 
variables  connected  by  M  links  and  calculate  the  upper  bounds  for 


Pr{  (j  A 
se 

there  is  onlv  one  such  set  x 


U  A  '  where  A.  is  as  defined  in  Theorem  A-TT.  In  the  above  example 
et  1 

M  M  M 


k.  x.,.  The  (A-Il)  upper 


bound  for 


Pr{ 


i  =  3, 5 ,6, 7 


(y )  ;•  is 


.  c , 


*  I"  [  Pr  -’A  j  ■  •  Pr  (A^  :  A^ }  •  Pr  {A^  A^  >  •  Pr  { Aj.  j  Ap  ] 


k  C  •  .  C  - 


by  repeated  applications  of  A-II. 


Note  the  order  that  (A-II)  was  applied  does  not  matter.  For  instance  if 
M  M  M 

we  took  x.  -w  -v  x-.  to  get 

6  o  3  7 


-  7t7\  C  ,CtP  OrP  C  >  C 

Pr;  (A.);  <  l-[  Pr  •.  A, }  •  Pr  lA^T  j  A.I  j  »Pr  (A  !  A,  ) *Pr {aV  A^} ] 

.  O  ,  ,  T  1  —  I  ->J  6  5’ 

1=3, 5, 6, 7 

the  reader  can  check  that  gives  the  same  value  as  *. 

So  let  i  be  the  (A-II)  upper  bound  for  Pr{  [  i  A.}. 

i-3.5,6.7  1 

The  next  step  is  to  use  Theorem  (A-I)  to  compute  upper  bounds 
using  those  calculated  in  the  first  step  and  incorporating  the  S 
links.  In  the  above  example  this  will  give  us 
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Pr{  u  (A  )} 
i=1.2 


<_  Pr{A^}  +  Pr{A9)  -  PriA^nA^)  by  theorem  (A-I) 


and 


Pr{  tj  (A  )} 

i=3 , 5 , 6 , 7 , 4 , 8 


<  v  +  Pr(A,  }+PrtAQ }-Pr{A, nA  }-PrlA0r>A  }  by  a 
—  4  o  4  1  o  / 

slight  modification  of  Theorem  (A-I)  . 


Let  [L  be  the  uoper  found  for  Pr{  U  (A.)}  which  in  this  case 
1  '  1=1,2  1 
was  Pr{Aj}+Pr{A9}-Pr  A.-A  }  and  let  ^  be  the  upper  bound  for  the 

Pr{  ii  (A.)}  which  was  derived  above. 

i-3,5,6,7,4,8  1 

The  last  step  is  to  incorporate  the  I  link  between  and 

by  using  corollary  (A-IIIC)  to  give  us  the  upper  bound  for 

8 

Pr{  .  ,  (A.)  >  <  1  -  (1-+  )  (1-ED2)  . 

i=l  1 

It  should  be  noted  that  for  real  cases  where  a  tree  has  been  determined  , 
assignment  of  S  and  I  links  is  arbitrary,  whereas  M  links  can  be 
assigned  only  if  the  joint  distribution  of  the  entire  group  connected 
by  M  links  meets  the  MTP0  conditions  of  theorem  (A-Il). 

It  should  also  be  mentioned  that  the  M  links  were  handled  first 
in  step  1  because  theorem  (A-II)  can  only  be  used  wich  M  links. 

There  is  no  way  to  incorporate  other  types  of  links  into  the  results. 

The  S  links  were  then  incorporated  before  the  I  links  were  since  doing 
so  gives  superior  results  (lower  upper  bounds). 
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Example  M-two. 

This  is  another,  different  approach  which  could  be  used  to 
combine  the  methods  together  to  obtain  an  upper  bound.  Let  us 
assume  x  is  a  4 x  1  vector  with  a  N(0,T)  distribution  and 
that  a  directed  tree  with  the  following  links  has  been  set  up 


X1  "  X2  "  X3  "  X4  • 


Let  Pr{A.}  have  the  usual  definition  for  i  =  and  assume 

1 

that  the  ( ; x  j , | X£ I , i x^ | )  are  jointly  MTP„  and  the  ( | x^ I , I | , I x 
are  jointly  MTP?  but  ( ! x^ j , | ' , ! x^ ! , ! x^ |)  is  not  jointly  MTP2. 
Then  by  Theorem  (A-II)  an  upper  bound  for  PrfA^-A^- A^'f  is 

-T  =  1-  [  Pr  {A^ }  *Pr  {A^’A^  }  /Pr  (a^  J  *  Pr  {a^'A^}  /  Pr  {A^ }  ]  * 


and  by  Theorem  (A-II)  an  upper  bound  for  Pr{A2-A^ jA^J  is 


i  =  l-IPriA^l-PrtA^'A^)  /Pr (A^ > -Pr  (A^Agl/Pr tA^} 


Using  and  ^  along  with  corollary  A-IC  and 


=  ,  S2  A 2 » A 3 , A ^ ,  S3  ^2*^3* 


^1  ^  ^2  ^1  ’^2*^3  *^4 


gives : 


Pr'  u  A.1-  <  +  ?2  '  Pr’A2  A3~ 

i=l 


Pr{A?.AC) 

Note  Pr{AC|AC}= - - — -1— 

1  3  Pr{AC} 
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It  would  be  nice  if  it  were  possible  to  find  the  best  tree  T 
and  method  combination  set  up  which  gave  the  lowest  upper  bounds. 
To  do  so  would  require  a  quadratic  program  which,  except  for  small 
n  would  be  computationally  impractical  to  solve. 


It  might  be  possible  to  develop  a  "good"  procedure  to  combine 

the  three  approaches  in  such  a  way  to  give  a  lower  upper  bound 

than  does  the  Dunn-Sidak  method  when  I  is  not  diagonal  and  hopefully 

give  an  upper  bound  which  is  close  to  the  lowest  possible  upper  bound 

achievable  from  any  procedure  combining  the  methods. 

The  author  believes  that  there  is  one  particular  situation  where 

n 

one  can  find  the  lowest  upper  bound  for  Pr{  ,  ,  A.)  nossible  from  com- 

i  =  l  1 

bining  the  three  procedures  and  show  that  this  method  is  indeed  the 
lowest  possible  upper  bound.  That  situation  is  when  the  absolute 
values  of  all  the  elements  in  x  (i.e.  \ ; , j x9 I , . . . , ] x^ j )  are 
jointly  MTPt 


Hypothesis  M-l.  Let  x  -  N'(0,~)  if  jx;  is  also  jointly  MTP  i.e. 

for  some  diagonal  matrix  D  wTith  elements  +1  such  that  the  off 

diagonal  elements  of  -DID  are  all  nonnegative,  then  the  upper  bound 
n  n 


for  Pr  {  A.}  or  Pr{  x.d(-c.,c.)}  given  bv  (A-II)  from  using 

^  =  ^1  .  1  11 

1  ~  ■*-  n 

any  tree  T  is  low^er  than  the  upper  bound  for  Pr{  ij  x.d(-c.,c.)} 

i=l  1  11 

given  from  using  any  combination  of  A-I,  A-IC  A-I,  A-1II  and 
A-IIIC  on  that  same  tree  T. 

The  proof  would  be  quite  complicated  since  there  are  so  many 
ways  to  combine  the  different  nethcds  together  which  must  be  considered. 
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Given  that  hypothesis  M-l  is  true,  then  to  get  the  lowest 

n 

possible  upper  bound  for  Pr{  \\  A.}  from  any  combination  of  the 

i=l  1 

three  approaches,  we  use  the  (A-II)  approach  entirely  on  T'  the 
best  tree  derived  in  remark  2  of  theorem  (A-Il)  ar.d  theorem  (AP-1) 
in  the  Appendix. 
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V.  EXAMPLES  OF  APPLYING  THE  COMBINED  METHODS. 


Example  One.  Consider  the  MA(5)  model.  Let  x  be  iia  N(0,1)  and 


let  y.  =  E  x .  ..  Then  P  ,  =  max( - L,0)  and  a  =  /5. 

J  i=l  J+i  yj,yj  5  y  j 


Let  the  set  R  =  '->’p 3 7 ’^8 ’ y9'^  an<^  suPPose  the  researcher 
wants  to  know  upper  bounds  for 


A) 

Pr{  'J 
jcR 

iv.!  , 

1.96 

/5} 

B) 

Pr{  U 

jeR 

iyj  > 

2.25 

/J} 

C) 

Pr{ 

jtR 

A 

*<“> 

>. 

2.50 

v'S-! 

D) 

A  lower 

bound 

for 

c  such  that  Pr{  1  ; 

y,  ;  >  c>  <  .10 

jcR 


Clearly  the  best  tree  T'  is  y ^-N^-y^-y y-Vg-yg  and 
?12  =  P9g  =  Pyg  =  P89  =  .80  while  P^y  =  .20.  It  can  be  shown  that 
for  any  continuous  group  of  size  3  or  more  on  T'  that  there  is  no  MTPy 
density  so  we  are  dealing  only  with  S  and  I  links.  Below  answers  to 


A,  B  and  C  are  shown  using  the  procedure  of  example  M-one  in  chapter  IV  with 


(i) 

all 

links 

as 

I  links 

(ii) 

all 

links 

as 

S  links 

(iii) 

>Tv 

2’y2~- 

3  ’ v 

_-y0,v  -v  as  S  links  y  -y_  as  an  I  link 
7  "  0  0  ‘  y  -it 

METHOD 

(i)  All  I  links  (ii)  All  S  links  (iii)  v^-Vy  1  link 

all  other  S  links 


A) 

Pr{  ;  |v 
j  -  R  '  J 

1  >1.96»^5> 

.2649 

.2051 

.2010 

B) 

Pr {  ;  |y . 

jcR  J 

!>2.25v/5) 

.1380 

.1076 

.1057 

C) 

P  r  {  i  [  v  . 

jcR"  j 

[  >2.50/5} 

.0722 

.0569 

.0564 
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The  combined  method  (iii)  works  better  than  does  either  the 
pure  1  method  or  the  pure  S  method  although  the  improvements  from 
(iii)  over  (ii)  are  not  spectacular.  The  upper  bound  ratio  of  method 
(iii)  over  method  (ii)  ranged  from  .991  for  C)  to  .980  for  A). 

The  answer  to  D  using  methods  (i),  (ii)  and  (iii)  respectively 
is 

(i)  2.383 

(ii)  2.283 

(iii)  2.275 


Example  Two. 


,xfe)  N (0, .  )  where 


1 

:  =  .9  l 

.9  .9  1 

.9  .9  .9  1 

.5  .9  .9  .9  1 

.3  .5  .9  .9  .9  1 


Let  the  researcher  be  interested  in  upper  bounds  for 

6 

A)  Pr {  x. '  "  1.96' 

i  =  l  1 

6 

B)  Fi-;  x .  ■  2.25  - 

.  i 


C)  Pr(  !x.  1  "  2.50^ 
i=l  1 


D)  A  lower  bound  for  c  such  that  Prt  :  x. ;  >  c 1  <  .05  . 

i=l  1 


Clearlv  a  best  tree  T’  is 
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X. -X-— x0-x . -xc-x , 

1  2  3  4  5  6 

It  turns  out  that  while  x  itself  is  not  jointly  MTP^  the  sub¬ 
vectors  ( 1  x^  I,' x^  ! )  and  (  ] x^  x^  x^  ',J  1  )  are  each  jointly  MT  P  0 . 

So  it  is  possible  to  use  an  algorithm  similar  to  that  of  example  M-Two 
in  chapter  IV  to  get  an  upper  bound.  We  shall  call  this  method  (*). 

(*) 

First  use  (A-II)  to  get  upper  bounds  for 

4  6 

*  =  Pr  {  ■  x .  '  >  o 1  and  •  =  Pr  •.  x.  >  c  ■  . 

1  .  ,  i  2  .  ,  l 

i=l  i=3 


Then  use  (A-i-C)  to  give 


6  4 

Pr {  x .  or  ■  •  +:  -  pr-  x. 

.  ,  l  -12  .  ,  i 

i=I  i=3 


Below  answers  to  A,  B  and  C  are  shown  using 


(i)  All  links  as  I  links 

(ii)  All  links  as  S  links 

(iii)  Method  (*) 

METHOD 


6 


A) 

Pr-  ■ 

i=l 

6 

X  . 

1 

>1.96 

B) 

Pr { 

i  =  l 

6 

x .  1 
1 

>2.25 

D) 

Pr' 

i  =  l 

1 

x  . 

1 

1  >2.50 

(i)  All  I  links  (ii) 
.  2643 

.1380 

.0722 


All  S  links 

(iii)  ■ 

.1518 

.1491 

.0792 

.0784 

.0423 

.0421 

2  3 


Again  the  improvement  of  method  (iii)  over  method  (ii)  was  not 
spectacular.  The  upper  bound  ratio  of  (iii)  over  (ii)  was  .995 

in  C ,  ,990  in  case  B  and  .982  in  case  A. 

The  answer  to  D  using  methods  (i) ,  (ii)  and  (iii)  respectively 
is 

(i)  2.635 

(ii)  2.436 

(iii)  2.433 


Example  Three.  This  is  based  on  problem  E  on  page  230  of  Drapper 
and  Smith.  A  linear  regression  is  de' eloped  to  predict  corn  yield 
per  acre  in  Iowa.  Five  of  the  coefficients  in  the  final  equation 
are  1.  Year,  2.  duly  Rain,  3.  Julv  Temperature,  4.  August  Rain 
and  5.  August  Temperature.  Ignoring  the  constant  coefficient  bO; 
call  the  five  coefficients,  bl,  b2,  b3,  b4,  b5,  and  the  estimated 
coefficients  bl,  b2,  b3,  b4  and  b5  respectivelv.  The  correlation 
matrix  of  b  is 

i 


The  best  tree 


-.1160 

1 

.4870 

.  344  5 

i 

-.1725 

.1215  - 

■.0914 

1 

-.2076 

.33o2  - 

■.3042 

.3339  1 

turns 

out  to  be 

'  1 

A 

—  ■.  b3 

b2 

/*  r’* 

b  c  -  •  b  . 

J  *4 

The  distribution  of  b  is  not  MTP0  but  the  distribution  of 

~  *- 

A  * 

b  9,  b^  and  b^  is  MTP?  however. 

Let  us  assume  that  the  researcher  wants  to  have  upper  bounds  for 
5 

A)  Pr;  ;  j  b.  i  b,  +  1.96  ' 


i=l 

5 


i  i 


B)  Pr'  ii  b .  4  b .  +  2.25 

i4  1  1  - 

5 

C)  Pr'  ii  b  .  4  b  .  +  2 . 50 

i  l  — 


b. 

l 


b. 

l 


b. 

l 


and 


D)  A  lower  bound  for  c  such  that  Pr;  b^  i  b^  +  ccg  :  .10 

i  =  l  1 


and  also  the  researcher  wants  a  lover  bound  for  c  so  that 
5  * 

Pr'  'u.  (b.~c"  )•  -  .10.  We  can  calculate  these  bounds  bv  using  the 

-=1  1  "  1  h.- 


procedure  ^.r.  example  M-one  section  IV  with  (i)  all  S  links,  (ii)  all  I 

•—o  I 


S/TN  1,-^.  M  M  ST' 

links  or  (iii)  the  following  links  (  b,; — ^  bQ  — !  bn  ' —  b_  — ^b, 


The  upper  bounds  for  A,  B  and  C  are 


A)  Pr-  b.tfb.il  ■' 

1  1  b 

i 


i  =1 


o  - 

B)  Pr;  '  b  .  ^  b  .  :  2 . 2  5 
.  ,  1  1 
i  =  l 


h. 

l 


i=l 


METHOD  (i) 
All  I  links 

.2262 


.1164 


5  - 

C)  Pr-  ■  b . ^b . i 2 . 50  c ,  >  .0605 

1  1  b 

i 


METHOD  (ii) 
All  S  links 

.  2245 


.1135 


.0589 


METHOD  (iii) 
.2156 

.1117 

.0586 
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The  upper  bounds  for  D  using  (i),  (ii)  and  (iii)  respectively  are 

(i)  2.319 

(ii)  2.310 

(iii)  2.306 

The  probability  ratio  of  method  (iii)  to  method  (ii)  was  .995  in  case 
C,  .984  in  Case  B  and  .960  in  case  A. 

Analysis . 

In  the  three  examples  of  section  V,  the  merged  methods  produced 
lower  upper  bounds  for  probability  than  did  any  of  the  single  methods 
used  alone.  In  all  three  examples  using  only  S  links  was  the  best 
competitor  being  from  .960  to  .995  as  efficient.  One  difficulty 
involved  with  merging  the  methods  is  that  it  is  not  simple  to  figure 
out  "how"  to  set  up  the  combined  links.  It  would  appear  difficult 
to  come  up  with  a  simple  universal  algorithm  to  do  this  meaning  that 
one  would  have  to  do  this  oneself  by  making  educated  guesses  as  to 
which  combination  of  methods  to  use. 


26 


VI.  TABLES  OF  RESULTS  FROM  COMBINED  (A-I)  AND  (A-III)  METHODS. 


At  this  point,  we  will  try  to  explore,  in  some  systematic  fashion, 
when  combining  the  methods  gives  improved  results  over  using  any  of  the 
methods  singly.  For  simplicity.  Method  (A-II)  the  conditional  multipli¬ 
cative  method  is  not  being  considered.  If  (A-II)  links  could  be  included 
in  the  comhined  methods  studied  in  the  next  few  tables  one  would  expect 
that  even  lower  upper  bounds  would  be  obtained.  The  methods  will  be 
combined  together  in  the  fashion  of  example  M-one  in  Chapter  IV. 

Tables  A,  B  and  C  look  at  various  combinations  of  S  and  I  links 


compared  to  the  best  of  using  all  S  links,  or  using  all  I  links,  for 

n 

calculating  upper  bounds  of  Pr;  x.^(-c,c)}  where  n= 4,6,10 


respectively, under  different  correlation  structures  for  the  links.  In 


all  three  tables  c  is  set  equal  to  1.96  and  2.50  while  the  best 

tree  T'  is  assumed  to  be  a  string  of  the  form  (x, -x„- • • • ~x  ). 

14  n 

SIS  SSI 

In  table  A,  link  patterns  of  the  form  x^x^x-^-x^  and  x.^x^x.j-x^ 

S  S  S  III 

are  compared  to  the  best  of  and  x.-x,-x,-x,.  The  combined 

1234  1234 

methods  do  give  better  results  than  does  the  pure  S  method  or  the  conser¬ 
vative  assumption  of  independence  v?hen  !P..  <■  .3  for  x.  linked  to  x. 

i.l  -  i  J 

by  the  I  link.  The  combined  methods  do  better  when  c  =  2.50  having 
upper  bound  ratio  values  as  low  as  .953  compared  to  pure  methods  than 
they  do  when  c  =  1.96  having  upper  bound  ratios  by  factors  as  low  as 


.977  compared  to  pure  methods. 

S  S  I  S  S 

In  table  B,  link  patterns  of  x^-x^-x^-x^-x^-x^  and 

SSSIS  SSSSS 

x,-x,-x,-x.-x,-x,  are  compared  to  the  best  of  x,-x„-x--x.-x.-x 
123456  12340 
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and  x  -x0^x  ix  ix, .  The  combined  methods  give  better  results  than 
12  3  4  5  6 

do  the  pure  methods  when  |p..[  <  .4  for  x,  linked  to  x.  bv  the  I 
link.  The  relative  upper  bound  ratios  of  the  combined  methods  to  the 
pure  methods  were  somewhat  lower  (better)  when  n  =  6  than  when  n=4.  When 
n  =  6  the  ratio  of  upper  bound  combined/best  uDper  bound  pure  methods  were 
as  low  as  .955  for  c  =  1.96  and  as  low  as  .985  for  c  =  2.50. 

In  table  C,  the  link  patterns  of  : 


, S  I  S  I  S  I  S  I  S„ 
(l)  xfx2-x3'!t4'x5"xrx7'5tg'x9'xio 

. .  .  s  SSiSISISS^ 
(n)  X2-X2“X3_X4_X4"X6"X7_X8-X9  10 


, . . . .  S  S  S  S  I  s  sv  sv  s 
(ill)  x1-x2-x3-x^  x5-x6-x7-xg  Xg  x10 


are  compared  to  the  best  of 


S  S  S  S  S  S  S  Sv 
xi“x2"x3-x4‘x5-x6_x7''x8  9  10 


and 


I  I 

xrx2_x 


i  i 
3~x4  x 


I  I  I  I 
5~X6_X7'X8 


9~X10 


The  combined  methods  give  better  results  than  do  the  pure  methods 

when  P  '  .8  for  c  =  1.96  and  |P  .  !  1  .6  for  c  =  2.50.  The 

ij  -  U 

probability  ratio  of  combined  method/best  pure  method  were  as  low  as 
.928  for  c  =  1 . 96  and  as  low  as  .977  for  c  =  2.50. 

Summarizing  tables  A,  B  and  C  the  following  tendencies  were  observed: 
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1.  The  combined  methods  were  better  with  respect  to  the  pure 
methods  as  n  became  larger.  It  is  felt  that  this  will  continue 
to  be  true  for  larger  values  of  n  than  were  looked  at  here. 

2.  For  |p..|  not  too  close  to  zero  where  x.  is  linked  to 

ij  i 

Xj  by  an  I  link  in  the  combined  method,  the  combined  methods 
became  better  with  respect  to  the  pure  methods  as  |p  |  became 
smaller . 

3.  The  combined  methods  were  better  with  respect  to  the  pure 
methods  for  the  smaller  value  of  c  =  1.96  versus  the  larger  value 
of  c  =  2.50.  It  is  felt  that  this  will  hold  in  general  for  all 
values  of  c. 

The  question  remains  as  to  whether  there  are  applications  where 
it  is  worth  the  extra  effort  to  compute  the  merged  results  to  get 
the  lower  upper  bounds  than  are  given  by  the  pure  methods. 
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TABLE  A 


UPPER  BOUNDS  FOR  Pr{  ft  x.£(-c,c);  OBTAINED  FOR  N=4, 

i  =  l  1 

WHEN  THE  BEST  TREE  IS  A  STRING  x^x^Xj-x^  BY  USING 
DIFFERENT  COMBINATIONS  OF  S  AND  I  LINKS 


ions 

along  the  various  links 

c 

Best  of 

=  1.96 

Smallest 

c  = 

P . . =1  Best  of 

S1^  all  S  and 

all  I 

2.50 

Smallest  P . . 

1~X2 

X2~X3 

X3_X4 

all  S  and 
all  I 

Rest  are 

Rest  are  S13 

.9 

.5 

.9 

.1314 

.1357 

.0355 

.0365 

.8 

#  5 

.8 

.1469 

.1501 

.0397 

.0405 

.9 

.4 

.9 

.1340 

.1357 

.0360 

.0365 

.8 

.4 

.8 

.1495 

.1501 

.0402 

.0405 

.6 

.4 

.6 

.1684 

.1675 

.0453 

.0451 

.9 

.3 

.9 

.1358 

.1357 

.0362 

.0365 

.8 

.3 

.8 

.1513 

.1501 

.0405 

.0405 

.  7 

.3 

.  7 

.1621 

.1600 

.0432 

.0432 

.6 

.3 

.  6 

.1703 

.1675 

.0451 

.0451 

.9 

.9 

_  2 

.1372 

.1362 

.0365 

.0365 

.9 

,  2 

.9 

.1372 

.1357 

.0365 

.0365 

.9 

.1 

.9 

.1379 

.1357 

.0366 

.0365 

.8 

.1 

.9 

.1455 

.1430 

.0387 

.0385 

.  5 

.5 

.5 

.1721 

.1731 

.0457 

.0464 

.5 

.4 

.5 

.1747 

.1731 

.0487 

.0464 

.4 

.  4 

.4 

.1799 

.1779 

.0487 

.0474 

.3 

.3 

.  3 

.1854 

.1811 

.0487 

.0480 

.3 

O 

.  3 

.1854 

.1811 

.0487 

.0480 

2 

.2 

2 

.1854 

.1837 

.0487 

.0484 

4 

Note  for  N=4  the  upper  bound  for  Prf  1  x 

i=l 

conservative  independence  assumption 

(i)  . 1834  for  c  =  1.96: 

(ii)  .0487  for  c  =  2.50: 

. i (-c , c) :  given  by 

is 

the 

30 


TABLE  B 


UPPER  BOUNDS  FOR  Pr{  (-c , c)  }  OBTAINED  FOR  N=6, 

WHEN  THE  BEST  TREE  IS  A  STRING  x,-x„-x~-x  -x„-x,  BY 

1  2  3  h  j  6 

USING  DIFFERENT  COMBINATIONS  OF  S  AND  I  LINKS 


Correlations  along  the  various  links  c  =  1.96  c  =  2.50 


1  X2 

x2-x3 

x^-x 

4  X4'X5 

VX6 

Best  of 

All  S  and 
all  I 

Smallest  P . ,=I 

Rest  are  S 

Best  of 

All  S  and 
all  I 

Smallest 
Rest  are 

.9 

.9 

.5 

.9 

.9 

.1723 

.1733 

.0475 

.0482 

.9 

.9 

.4 

.9 

.9 

.1749 

.1733 

.0479 

.0482 

.8 

.8 

.4 

.8 

.8 

.2057 

.2011 

.0562 

.0563 

.9 

.7 

.4 

.  7 

.9 

.2010 

.1970 

.0549 

.0554 

.6 

.5 

.  4 

.5 

.6 

.2499 

.2401 

.0669 

.0666 

.4 

.4 

.4 

.4 

.4 

.2665 

.2545 

.0702 

.0698 

.8 

.7 

.  3 

.7 

.8 

.2183 

.2107 

.0594 

.0590 

.8 

.8 

.8 

.3 

.8 

.2075 

.2019 

.0566 

.  0564 

.3 

.3 

.3 

.3 

.3 

.2649 

.2607 

.0722 

.0711 

.3 

.3 

.1 

.3 

.3 

.2649 

.2607 

.0722 

.0711 

Note 

for  N=6 

the  upper 

bound  for  Pr 

6 

{  )  x .  i  ( -  c ,  c ) } 

given  by  the 

i  =  l 

conservative  independence  assumption  is 


(i) 

.2649 

for 

c  =  1.96 

(ii) 

.0722 

for 

c  =  2.50 
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UPPER  BOUNDS  FOR  I’rl  i  .j  x.i/(-c,c)}  OBTAINED  FDR  N'=6  WHEN  THE  BEST  TREE  T  IS  A  STRING 


CD 

-X 


u 

z 

M 

t—J 

CO 

X 

1  X 

>- 

oo  xc 

02 

X  c 

►  *H 

r—» 

r-«.  *-h 

a. 

X 

1  co 

W 

vC 

C3 

X  -UJ 

»  X 

< 

tn  <y 

> 

X  0 

1 

:__ 

1 

<r 

Z 

X 

< 

» — 

ro 

CO 

X 

z 

1 

Q 

CNJ 

CJ 

X 

< 

COHnNfnNMCMT' 

lOvjmrNWNnrsO 

O^invCrNCOO'C'O 

nnnnnfonrir) 


r^^rr-OOcNr^oor^ 
oooO'JOncs1  <rmvor^ 
OOCTnOO^H^^H’-H'H 
O  O  r-H  ’“I  «— I  t— I  »— I  »— !  »— t 


X 

oc  -X 
X  c 


II 

u 


X  4— < 

!  X 

lT'  O 

X  OL 


'sCOOC''1 — I  H  IT  vj  M 
inHinrsivDrsNjCH 
OCCM^vOr^oOOOO 


X 


I 

CN 

X 


X 

-X 


X 

.X 


x  ^  h  vC  ^  n  n  rj  n 
inXsrcc^-J'J'' 
^ro'c^v.cccCT'OCC 


C  *—  LT|  c>  lT>  CNJ  CM  CM  CJ  CN 

CO  w  C'*"'  CNJ  *-H  00  *”H  *~H  »-H  *-H 

•u  mo^rOOCOCC 

x  i—  T3  cNjr^r^m<r<r^r<r<r 

Q  C2  ••••••••• 

^  tr  c 


^  CO  P'  r  LTi  <J  Csl  H 


c 

u->  O'vOCNjinCNJ'^O'lTi 
CM  cc  O'  C  c  —•  — ^  * — >  ’ — 1  ■ — 1 

ii  *""• . 

u 


<T  O'  ^  LO  CC  UN  LO 
r— i  r^r^.LOO^vCr^.r-- 

^  ^  c  - - I  - - 1 

'  C  *”4  l—1  —* *  *“ -i  r-H 


C"»  I— ‘  lH  o  <r  iT>  i/N 

vOit.  r>-vCr4vCNiNr^ 
nC  CO  C  -  C  *”“1  ’“4  1  ”4  <— • 
O  O  O  «-h  ^  » — 1  *““•  *“H  •—< 


c.  S'  cc  n  c  to  c  ^  p;  *4 


32 


VII.  TERMINOLOGY 


x  -  is  an  nxl  vector  of  random  variables  with  a  N(0,v)  distribution 

P..  -  is  the  correlation  between  x.  and  x. 
iJ  1  1 

-  is  the  event  x^  is  not  in  the  interval  (-c^,c^) 

S.  -  is  a  set  of  integers  representing  corresponding  numbered  events 
from  A^,A2»...,A  ^  or  equivalently  corresponding  numbered 

variables  from  x^,X2,...,x^  ^ 

(mA  n)  -  is  the  minimum  of  m  and  n 

Pr{ • }  -  is  the  probability  that  {•}  occurs 

T  -  is  a  tree  connecting  the  events  A^,A2,...,An  or  equivalently 

the  variables  x. ,x~ , . . . ,x 
12’  n 

e..  -  is  a  directed  link  from  x.  to  x.  or  equivalentlv  from 
ij  l  J 

e .  . 

A.  to  A,  e.g.  x  22.  x 
1  J  i  J 


S  link  -  an  S  link  from  x.  to  x.  indicates  that  the  upper  bound 

i  J 

operation  which  will  be  performed  between  these  two  variables 
(events)  is  of  the  intersection  subtraction  (Theorem  A-I)  type. 

M  link  -  an  M  link  from  x.  to  x.  indicates  that  the  upper  bound 

i  j 

operation  which  will  be  performed  between  these  two  variables 
(events)  is  of  the  conditional  multiplicative  (Theorem  A-Il)  type. 


I  link  -  an  I  link  from  to  x  indicates  that  the  upper  bound 

operation  which  will  be  performed  between  these  two  variables 
(events)  is  of  the  conservative  independent  subunit  (Theorem  A-III) 


s  -  an  upper  bound  to  Pr;.  A.'. 

i  =  l  1 

A  B  -  Pr  A  B.  means  the  probability  of  A  given  B. 
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VIII.  APPENDIX. 


Below  is  presented  the  theorem  to  find  the  "best"  tree  with 
which  to  apply  method  (A-I)  or  method  (A-II)  for  .he  case  m=2. 

Note  that  a  maximal  tree  is  one  which  has  the  maximal  number  of  links 
connecting  n  elements  while  forming  no  loops,  (n-1  links). 


Theorem  AP-1.  Let  x  'u  N(0,E)  and  let  P..  be  the  correlation 
-  tj 

coefficient  between  x.  and  x..  Let  T  be  any  maximal  tree 

1  1 

connecting  the  individual  components  of  x.  Let  0^.  be  the  order 
statistics  of  jP_|  for  e_tT.  Set  up  a  tree  T*  by  performing 


the  following  process 

which  maximize  !p..; 

ij 


n-1  times:  Link  x.  to  x,  for  the  i,i 

i  J 

among  all  unlinked  i,j  whose  linking  does 


not  create  a  circuit.  Then  this  tree  will  have  the  following  property 


(0T,)1  >_  f°r  al-l  other  trees  T  and  all  j  =  l,...,n. 

Proof .  The  reader  is  referred  to  Hoover  (1986)  or  Kruskal  (1956). 


Remark:  Since  Sidak's  theorem  gives  that  Pr(A^'A^)  is  a  nondecreasing 

function  of  'p..!  the  above  theorem  implies  that  T'  is  the  best  tree 
ij 

on  which  to  apply  methods  (A-I)  and  (A-II).  Note  than  when  using  a 
combined  method  for  any  given  situation,  T  may  not  be  the  best  tree, 
but  should  still  be  a  good  tree  to  use. 
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